Due to the shortcomings of large amount of computational grids and slow convergence rate in the numerical simulation of the complex flow, the parallel algorithm of multiple Cartesian grid generation based on three-dimensional geometry was proposed for Lattice Boltzmann Method (LBM). And multi-grid parallel LBM was developed based on the parallel algorithm of grid generation. The proposed algorithms can efficiently reduce the amount of computational grids and improve the convergence by coupling with the grids with different sizes. Numerical results also show that the proposed parallel algorithms have considerable scalability.